A mechanism for unipolar resistance switching in oxide non-volatile memory devices 
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Building on a recently introduced model for non-volatile resistive switching, we propose a mecha- 
nism for unipolar resistance switching in metal-insulator-metal sandwich structures. The commuta- 
tion from the high to low resistance state and back can be achieved with successive voltage sweeps of 
the same polarity. Electronic correlation effects at the metal-insulator interface are found to play a 
key role to produce a resistive commutation effect in qualitative agreement with recent experimental 
reports on binary transition metal oxide based sandwich structures. 

PACS numbers: 85.30.Tv, 85.30.De, 73.40.-c 



Resistive switching phenomena controlled by external 
voltages on metal-insulator-metal (MIM) sandwich struc- 
tures is currently attracting a great deal of attention 
due the potential applications for non-volatile memory 
devices. The switching is an abrupt change between a 
highly resistive state (off state) and a conductive state 
(on state) driven by the applied voltage. 

In recent years several MIM structures where the di- 
electric is a transition metal oxide (TMO) have been re- 
ported 0,0,0]. After electro-forming, i.e., the applica- 
tion of large voltages in order to produce a permanent 
change in its electric properties, those systems exhibit 
large hysteresis in the I-V curves. In most cases reported 
so far the nature of the switching is bipolar, that is, if a 
system is switched from one state to the other by a pulse 
of a given polarity, then a pulse of opposite polarity has 
to be applied to switch it backP, 

Surprisingly, several recent experiments on simple bi- 
nary oxides (Ti02, NiO, CoO, etc.) demonstrated a novel 
switching effect d, 0, S H S S]- these systems the 
switching can be achieved by successive application of 
voltage pulses of the same polarity. While some qualita- 
tive ideas have been proposed, the origin of this so called 
unipolar switching phenomenon remains largely not un- 
derstood. 

Recently, we have proposed a phcnomcnological model 
that explains several aspects of the bipolar switching phe- 
nomena 0, [l3| • The goal of this letter is to build up on 
that basic model and identify what additional hypothesis 
are required to account for the unipolar switching effect 
within the same framework. We anticipate that the key 
requirement is that the current injection during the high 
to low resistance transition (called SET) must be of short 
duration (see below). 

The phenomenological model introduced in RefsS, [13 
consists of an inhomogeneous dielectric medium con- 
nected to two electrodes (called top and bottom). The 



dielectric is conformed by two sets of small domains, each 
one in close physical proximity to the electrode interfaces, 
and also an additional large central domain that acts as 
a reservoir Carriers move from domain to domain 
under the action of an electric potential. The hopping 
path of carriers can be associated to filamentary con- 
duction across the dielectric. The domains may actually 
correspond to grains, dislocations, concentrated chemi- 
cal inhomogeneities, such as vacancies or other type of 
defects 0, [ul [13, 
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In the model 



they are simply characterized by the number of carri- 
ers that they can contain and their occupation at a given 
time. The model is further defined by the parameters 
that control the probabilities for the carriers to hop from 
domain to domain or to the electrodes. An important 
ingredient of the model is the effect of correlations [lo| . 
Since TMO are usually prone to strong correlation ef- 
fects ^1^, we assume that under certain conditions that 
can be tuned by the electric pulsing, the interface (small) 
domains might undergo a metal-insulator transition. As 
a result, the density of states of those domains open a 
charge (Mott) gap, and their associated hopping proba- 
bilities get greatly modified. Since the affected domains 
are located at the interface, this effect has a significant 
impact on the conductance of the whole system. 



Similarly as in RefilOl. here we adopt a simplified ver- 
sion of the model where the, in principle many, small 
"top" and "bottom" domains are represented by a single 
typical domain. Thus the model simply has only three 
domains: bottom, central and top. Since the systems 
are usually symmetric, following Refs. we adopt 

the same hopping rates, r*^^*, for the two interfaces [bot- 
tom(top) electrode - bottom(top) domain] and the same 
yint g^jj i-jjg internal hopping (domain - domain) . The 
model is thus defined by the following three coupled non- 
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linear differential equations [lo| . 



^ = r"*/(^)A^6«'(l -n') -r™*/(F)n^7Vt(l -n*),(2) 

^ = r«V(F)iv,n^(i-n*)-r-*/(yy^ (3) 

A^o! is the total number of states in electrode or do- 
main a, and n" is the occupation, a = te,he,t,c,h de- 
note "top-electrode", "bottom-electrode", "top", "cen- 
tral" and "bottom". The voltage dependent function 
fiy) actually gives a time dependence to the coefficients 
of the differential equation via the applied voltage pro- 
tocol V{t). This feature renders the system of equations 
very non-linear. A priori, the argument V of these volt- 
age dependent contribution to the hopping probabilities 
should be a fraction of the externally applied voltage, 
and a different one for each term in the equations above. 
However, for simplicity and to avoid a proliferation of 
model parameters, we shall assume that this dependence 
can be absorbed in the definition of the F's. Explicitly, we 
adopt f{V) ~ sinh(fcy) [3] where k depends on various 
material parameters, and again, for the sake of simplic- 
ity, we set A; = 1. Positive current is defined as carriers 
entering the bottom domain (or leaving the top domain) . 
As in RefH^ we assume that the electronic state of the 
domains may undergo a metal-insulator transition. In 
practice, this may occur due to injection (leakage) of par- 
ticles into (from) small domains, so that the occupation 
can trigger a density driven Mott metal-insulator tran- 
sition, which is common in transition metal oxides [l8| . 
Alternatively, it may also occur due to a local electro- 
chemical oxidation [71, 



]J, lia l20|, l2l| that is produced 



by the strong current densities during the electric puls- 
ing. The oxidation is likely to occur near the electrodes 
at the small domains J,, 22^ , and it would locally restore 
the originally insulating properties of the oxide. The phe- 
nomenon of local oxidation induced by a strong current 
has been observed in Ti films . 

Upon becoming insulating, the small domains open a 
charge gap A in their excitation spectrum. We consider 
this effect explicitly in the model by modifying the ap- 
propriate fiV) in Eqs.(l)-(3) to f.mtiV) = e-^/^f{V), 
where T is temperature (more details are given below). 
Following [lo| we assume that domains are insulating 
when 1/2-S < < 1/2+d, ie, when their occupation 
is near half-filling with about one electron per site. We 
set S = 0.1. 

Within this model, a mechanism that produces the 
unipolar switching is schematically shown in the sequence 
of panels presented in Fig.[T] Similarly as in the case de- 
scribed in Reffiol. the switching effect is mostly controlled 
by either one of the two small domains 0, 24 1 , which act 
as "electric faucets" Isl . Here we shall assume that the 
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FIG. 1: (color on-line) Schematic steps of the unipolar resis- 
tive switching process. The number of the panels corresponds 
to the numbers in the Fig[2] 



bottom one is the relevant one, while the other remains 
conductive all the time. 

Panel 1 schematically indicates the initial "electro- 
forming" step that is experimentally required. This pro- 
cess of soft breakdown is not described by our model. It 
involves the metalization of the initially dielectric bulk 
(central domain). We shall simply assume that form- 
ing sets the system in the state of panel 2, which is the 
initial state of our calculation. Thus, in panel 2, the sys- 
tem is in the conducting ON state (open faucet) where 
the small bottom domain (BD) is metallic with an occu- 
pation above 0.6. Upon application of a voltage pulse, 
the carrier hopping is enhanced and due to our assump- 
tion of F'^^* << r*"*, there is leakage from BD to the 
central bulk (panel 3). The leakage decreases the occu- 
pation of the BD until it reaches the threshold for the 
metal- insulator transition. At that point, the BD opens 
a gap in its band-structure and the hopping probabilities 
get reduced by a factor e~^l^ that enters in the renor- 
malization /(U) fmit{V), as described before. This 
metal-insulator transition is thus responsible of the RE- 
SET of the system that is now in the OFF state (panel 

4) . 

To switch the system back to the low resistance state 
(SET transition) a new voltage pulse is initiated (panel 

5) . At low applied voltage a very small current circulates 
since the BD is insulating (closed faucet). As V increases, 
a strong electric field develops across the BD. When the 
voltage at the electrode reaches a threshold value Ket, 
such that the Fermi level of the electrode aligns with the 
bottom of the empty conduction band of the BD, a sud- 
den injection of carriers occurs (panel 5). As this thresh- 
old is reached, the exponential factor in fmit(V) in the 
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FIG. 2; (color on-line) Top panel: Bottom domain occupa- 
tion as a function of the time. The numbers sequence refers 
to the panels of Fig. [T] The ON and OFF states are also 
indicated. Middle panel: Voltage protocol as a function of 
time. The SET and RESET voltages, T4et and Keset, are 
also shown. Bottom panel: Current entering the BD as a 
function of time. In this scale only the ON state current 
can be observed. The OFF current occurs between tr and ts 
cannot be distinguished from the x-axis. For clarity we show 
with dashed lines the RESET and SET transition times tr and 
ts for the third and fourth voltage pulses respectively. The 
model parameters are Nt = Nt, = 10^ = 10^ r'"* = 10"*, 



=10" 



and A/T = 10. 



first term of the rhs of Eq.([T]) disappears, and this V de- 
pendent coefficient goes back to f{V). It is important to 
realize that the second term of that equation remains un- 
modified as schematically shown in panel 5. At this point, 
the BD begins to increase its occupation and eventually 
the BD returns to the metallic state. This SET transi- 
tion (open faucet) is very fast; it can be associated to a 
new local electric breakdown at the BD, as it occurs at a 
relatively high applied V. Right after the transition, the 
large applied voltage produces a large current across the 
be/b interface and, therefore, a rapid increase of the oc- 
cupation of the small domain. Thus, to prevent the BD 
from leaking the carriers down to the central domain, the 
duration of the applied voltage pulse should be small. It 
is interesting to notice that this is reminiscent of the ex- 
perimental situation, where the SET transition has to be 
performed with a current control [2^ , which amounts 
to a rapid decrease of the applied V as the device be- 
comes conducting. After the BD returns to the metallic 
state (panel 6), the factor fmit{V) of the 2"'' term of Eq[T] 
also returns to f{V) (open faucet) and the system goes 



back to the ON state of panel 2; a new cycle may then 
begin. 

The SET process described in detail in the previous 
paragraph in which the system switches back from the 
OFF state to the ON state by appHcation of a short volt- 
age pulse of the same polarity as the initial reset pulse, is 
the key difference between this work and that of RefllO. 

We now substantiate the previous qualitative descrip- 
tion by solving the model equations following the consid- 
erations just described. In Fig[2]we show the occupation 
of the bottom domain as a function of time obtained 
by the action of six successive voltage pulses V{t), all 
with positive polarity. In the top panel of Fig [2] the con- 
secutive discharges and charges of the bottom domain 
are clearly observed. They can be correlated with the 
V{t) shown in the middle panel, where the thresholds for 
SET and RESET are indicated. The current I{t) enter- 
ing the BD is shown in the lower panel, where the high 
current injection after the SET transition is clearly ob- 
served. On the other hand, due to the linear scale that 
we use, the low current level between tr and ts cannot 
be distinguished. To correlate the successive instants of 
the calculation with the previous qualitative discussion, 
we include in Figl2] a number sequence that refers to the 
panels of FigH] We note that in the OFF state the oc- 
cupation of bottom domain is smaller than 0.6, but only 
by an exponentially small amount which cannot be dis- 
tinguished in the linear scale used. There is a certain 
freedom in the choice of the initial value of the occupa- 
tion of the BD [the value Ubit = 0)=0.62 in Figl2]. This 
arbitrary initial occupation would correspond to the ini- 
tial state after the electro-forming process which is not 
described by the present model. Nevertheless, the impor- 
tant observation is that after the first sequence of RESET 
and SET pulses occurs, the initial condition is lost and 
the occupation rib, along with the memory operation, be- 
come cyclic. 

To conclude, the behavior of our model is found to 
be in good qualitative agreement with the experimental 
reports of unipolar switching in transition metal oxide 
memories. We hope that our results provide a useful 
guidance for further experiments needed to improve our 
understanding of the basic switching mechanism in ox- 
ide non-volatile memories. The parameters of our model 
can be related to technologically important features such 
as switching speed and non-volatility. However their ex- 
plicit calculation using ab initio methods, which would 
open the way for a multi-scale modeling of memory de- 
vices, remains a great challenge ahead. 
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